Convergence to global consensus in opinion dynamics under a nonlinear voter model 
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We propose a nonlinear voter model to study the emergence of global consensus in opinion dynamics. In our 
model, agent i agrees with one of binary opinions with the probability that is a power function of the number of 
agents holding this opinion among agent i and its nearest neighbors, where an adjustable parameter a controls 
the effect of herd behavior on consensus. We find that there exists an optimal value of a leading to the fastest 
consensus for lattices, random graphs, small-world networks and scale-free networks. Qualitative insights are 
obtained by examining the spatiotemporal evolution of the opinion clusters. 

PACS numbers: 89.75.Hc, 87.23.Ge 



I. INTRODUCTION 

Mutual agreement, or consensus, is a fundamental phe- 
nomenon in social and natural systems. The dynamics of 
opinion sharing and competing and the emergence of consen- 
sus have become an active topic of recent research in statisti- 
cal and nonlinear physics fljj] . For example, a number of mod- 
els have been proposed to address how consensus can result 
from the evolution of two competing opinions in a popula- 
tion, which include the voter model ISyl], the majority rule 
model JUm, the bounded-confidence model [6] and the so- 
cial impact model [7]. Due to the high relevance of complex 
networks Hfl to social and natural systems, opinion dynamics 
have also been incorporated on networks l9l- l36ll such as reg- 
ular lattices, random graphs ll42ll . small-world networks ll43ll 
and scale-free networks fl44tl . 

Previous works have revealed phase transitions in opinion 
dynamics |1 [13, El El El 

and the emergence of global con- 
sensus H H where all agents share the same opinion. It has 
also been found that both the network structures 111 211 and the 
opinion updating strategies iflil \T]l [l8ll can affect the time 
for reaching the final consensus. InRef. iH, Yang et al. 
combined the majority rule model with probability p with the 
voter model with probability 1 — p and then measure the re- 
sulting consensus times on scale-free networks. They found 
that that the optimized ratio to minimize consensus time is 
around p = 0.72. In Refs. (Hffll, Yang et al. studied the ef- 
fects of heterogeneous influence of individuals on the global 
consensus. Each individual is assigned a weight that is pro- 
portional to the power of its degree, where the power expo- 
nent a is an adjustable parameter. Interestingly, it is found 
that there exists an optimal value of a leading to the shortest 
consensus time for scale-free networks, random networks and 
small-world networks. 

In the voter model, an agents follows the opinion of a ran- 
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domly selected neighboring agent. In the majority rule model, 
an agent follows the local majority opinion. For the voter 
model and the majority rule model, an agents absolutely fol- 
lows one opinion. However, in reality, there can be situations 
where an agent chooses one opinion with a stochastic proba- 
bility due to bounded rationality. In this Letter, we propose a 
nonlinear voter model in which an agent i selects one of binary 
opinions with the probability that is proportional to the power 
of the number of agents carrying this opinion among agent i 
and its nearest neighbors. The power exponent a is introduced 
to control selective probability. Our main finding is that, there 
exists an optimal value of a for which the convergent time 
for consensus is minimum. This phenomenon indicates that, 
a suitable preference of the local majority opinion (not fol- 
lowing it absolutely) can greatly accelerate the convergence 
process towards the final consensus. 



H. MODEL 

Our model is described as follows. For a given network of 
any topology, each node represents an agent. Initially the two 
opinions denoted by the values ±1 are randomly assigned to 
agents with equal probability. At each time step, agents syn- 
chronously update 114 ill their opinions according to the follow- 
ing rule. Agent i will select the opinion +1 with the probabil- 
ity 



where n + and n_ are the number of agents holding the 
opinion +1 and —1 among agent i and its nearest neigh- 
bors respectively, and a is an adjustable parameter. Sim- 
ilarly, agent i selects the opinion —1 with the probability 
p_ = n°:/(n°l_ + n"), where p_ = 1 — p + . 

The parameter a characterizes the degree of herd effect on 
consensus. For a > (< 0), agent i has a larger probability 
to select the local majority (minority) opinion. In the case 
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FIG. 1: (Color online.) Order parameter r\ as a function of t for 
different values of a on a 50 x 50 square lattice with periodic bound- 
ary conditions. Each data point is obtained by averaging over 2000 
different realizations. 



FIG. 2: Time to achieve consensus, T c , as a function of a for differ- 
ent lattice size N. The inset shows the optimal value of a, a opt , as 
a function of TV. Each data point is obtained by averaging over 2000 
different realizations. 



described as: 



of a = 0, agent i randomly selects the opinion +1 or — 1. 
For a = 1 , our model coincides with the voter model. For 
a — > oo, agent i absolutely takes the local majority opinion 
and our model is reduced to the model proposed in Ref. Il5ll . 



III. RESULTS AND ANALYSIS 

We first consider a square lattice with periodic boundary 
conditions. Following Ref. ifioll . we define an order parameter 
r\ as 



1 N 



(2) 



where <7j is the value of node i's opinion (+1 or —1) and N 
is the total number of nodes of the network. In general, we 
have < r\ < 1, and a large value of r\ indicates that one 
opinion in the system dominates the other. If the two opinions 
are equally probably, we have r\ — 0. When global consensus 
is achieved so that there is only one opinion in the system, we 
have i) = 1. 

Figure Q] shows the order parameter 77 as a function of the 
time step t for different values of a. For a=l, 1.1 and 2, r\ 
increases from to 1 with t. For a = 2.0, initially 77 increases 
faster than the cases of a = 1.0 and a = 1.1. Nevertheless, 
the system converges to the global consensus state faster for 
a = 1.1 than for a = 2.0. Another feature in Fig.Q]is that, 
for a < 1, the time required to converge to global consensus 
can be extremely long. 

Why does the system become so difficult to reach consen- 
sus for a < 1? According to the mean-field approximation, 
opinion dynamics in the preferential-selection model can be 



df>± = P+P- 
dt ~ p% + p c _ 

#(1-P+) 



_ P1P+ 
~ P+ + P- 



# + (i-p+)« p% + (i- P+ ) a 



p+i 1 - p+) 
% + (i-p+) 



-[pl-'-ii-p+T- 1 ], 



(3) 



where where p + and p_ are the fractions of agents holding 
the opinion +1 and —1 in the total population, respectively. 
In the case of a < 1, for p + > 0.5 (< 0.5), dp+/dt < 
(> 0). This indicates a negative feedback mechanism that 
prohibits the majority opinion from spreading over the whole 
population. 

We define the consensus time T c as the time steps required 
to reach the final consensus. Figure [2] shows that T c as a func- 
tion of a for different network size N. It can be seen that there 
exists an optimal value of a, hereafter denoted by a opt , which 
results in the shortest consensus time T c . We have found that 
a op t decreases as the network size N increases, as shown in 
the inset of Fig. [2] 

To understand the process of convergence to consensus, we 
study the evolution of opinion clusters. A opinion cluster is 
a connected component (subgraph) fully occupied by nodes 
holding the same opinion. Figure[3]shows the number of opin- 
ion clusters N c i and the normalized size of the largest clus- 
ter Si = S max /N as a function of t for different values of 
a, where S max is the size of the largest cluster. We see that 
eventually 7V C ; decreases to 1 and Si increases to 1 for large 
t. For small values of a, agents do not select the local major- 
ity opinion with large probability, thus it becomes difficult for 
them to form large opinion clusters. For with a = 1.0, at the 
beginning, N c i decreases and Si increases much more slowly. 
For a = 2.0, initially N c i decreases and Si increases more 
quickly than the cases of a = 1.0 and a — 1.1. However, 



3 





10° 10 1 10 2 10 3 10 4 10 5 

t 



Co" 





s 

/ ' / 
t, ' / 
/ / 
/ / 
•/ jT 




s 
/ S 


a=1.0 




y 

■ AT 


o=1.1 






a=2.0 


(b) 



I — — — — — I 

10° 10 1 10 2 1 3 1 4 1 5 



t 



FIG. 3: (Color online.) (a) The number of opinion clusters N c i and 
(b) the normalized size of the largest cluster Si as a function of time 
t for different values of a on a 50 x 50 square lattice. Each data point 
is obtained by averaging over 2000 different realizations. 



FIG. 5: Snapshots of opinion patterns on a 50 x 50 square lat- 
tice. Initially (t = 0), we set a subregion of 30 x 30 at the center 
where all agents have the opinion +1 (white), but the other nodes 
carry the opinion —1 (black). Initially the fraction of the opinion 

+ 1 is p+(0) = 9/25 = 0.36. (a) a = 0.95, p+(200) = 0.41, 
p+(500) = 0.4576, p+(700) = 0.474; (b) a = 1.1, p+(200) = 
0.2712, p+(500) = 0.1276, p+(700) = 0.0278; (c) a = 5, 
p+(200) = 0.3232, p+(500) = 0.2428, p+(700) = 0.1916. 




FIG. 4: Illustration of the configuration for the coexistence of two 
opinions for a — > 00. 



when there are only two or three opinion clusters remained 
in the system, N c i decreases and Si increases very slowly 
for a = 2.0, indicating that the merging of different clus- 
ters becomes difficult for large values of a. In particular, for 
a — » 00, the system cannot reach consensus in some cases. 
As shown in Fig. [4] two different opinions can coexist forever 
when a —5- 00. For moderate values of a (a = 1.1), large 
opinion clusters can be formed more rapidly than for low val- 



ues of a (a = 1.0), and the their emergence is facilitated than 
the cases of larger values of a (e.g., a — 2.0) as well. Thus 
the convergence time becomes minimum for some moderate 
value of a. It can also be noted that the evolution of Si for 
different values of a is similar to that of -q (see Fig.Q]). 

To assess how different values of a affect the evolution of 
the opinion clusters, we design a numerical procedure to in- 
vestigate the evolution of spatial patterns of opinion. Specifi- 
cally, in the 50 x 50 square lattice, initially we set a subregion 
of 30 x 30 at the center where all agents have the opinion +1, 
but the other nodes carry the opinion —1, as shown in Fig. [5] 
We see that for a < 1, e.g., a — 0.99, the boundary be- 
tween the two opinion clusters gradually become blurred and 
they tend to be well-mixed with approximately equal densi- 
ties. In this case, global consensus is less likely due to the 
difficulty to form opinion clusters. For moderate values of 
a, e.g., a — 1.1, the boundaries become irregular with time 
but finally the +1 cluster vanishes. For high values of a, 
e.g., a = 5.0, the boundaries are clear and the inside clus- 
ter shrinks but at a speed lower than that for a = 1.1. 

The existence of an optimal value of a for which global 
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FIG. 6: (Color online.) Convergent time T c as a function of a for ER 
random graphs, NW small-world networks and BA scale-free net- 
works. The average connectivities of ER, NW and BA networks are 
4. The inset shows the value of a op t as a function of the average 
connectivity (k) for the three types of networks. The network size 
is N = 3000. Each data point is obtained by averaging over 100 
different network realizations with 20 runs for each realization. 



consensus can be achieved rapidly is not restricted to the 
square-lattice structure. In fact, we have observed a sim- 
ilar behavior for complex networks including Erdos-Renyi 
random graphs (ER) 14211 . Newman-Watts small- world net- 
works JNW) 14311 and Barabasi-Albert scale-free networks 
(BA) 14411 . We also observe that, the optimal value of a tends 
to increase with the average connectivity (k) (see the inset of 
Fig. [6]). From Fig. [6] one can see that the topological structure 
affects the convergent time T c . Among the three complex net- 



works, T c for scale-free networks is shortest when (k) and a 
is fixed. 



IV. CONCLUSION 

In summary, we have proposed a nonlinear voter model 
to study the convergence to global consensus on lattices and 
complex networks. An agent i selects one of binary opin- 
ions with the probability that is proportional to the power of 
the number of agents carrying this opinion among agent i and 
its nearest neighbors. The power exponent a is introduced to 
govern selective probability. It is found that there exists an op- 
timal value of a leading to the shortest convergent time. We 
have explained such phenomenon in terms of the evolution of 
the opinion clusters. For too small values of a, the forma- 
tion of big opinion clusters is slow. For very high values of 
a, the merging of different opinion clusters becomes difficult. 
Taken together, the shortest convergent time can be realized at 
moderate values of a. Our results indicate that, the consensus 
would be quickly reached if agents suitably follow the local 
majority opinion. 
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